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I.  INTRODUCTION 

At  present  a  great  deal  of  Interest  is  being  given  to  the  practical 
control  of  large  space  structures  such  as  space  transportation  systems 
and  large  communication  satellites  [1].  The  control  task  is  normally 
thought  of  in  terms  of  maintaining  specified  shape  configurations,  orientation 
and  alignment,  vibration  suppression  and  pointing  accuracy,  to  name  a 
few.  Because  of  the  inherent  flexibility  associated  with  these  systems 
they  are  generally  analyzed  as  distributed  parameter  systems  which  creates 
difficulties  in  the  design  and  analysis  of  controllers  for  them.  Further, 
control  laws  derived  using  distributed  parameter  modeling  cannot  be  implemented 
efficiently  with  current  sensor  and  actuator  technology.  Modal  control 
techniques  have  been  developed  specifically  to  bypass  the  problems  associated 
with  distributed  parameter  theory  [2,3].  The  concept  of  modal  control 
is  built  upon  the  notion  that  certain  specified  system  modes  can  be  controlled 
by  appropriate  design  of  the  associated  closed-loop  eigenvalues.  This 
reduces  the  number  of  sensors  and  actuators  needed  to  effect  the  control 
of  the  structure.  However,  an  undesirable  phenomenon  referred  to  as  observation 
and  control  spillover  can  occur  if  the  number  of  sensors  and  actuators 
used  is  small.  Spillover  refers  to  the  phenomenon  in  which  energy  intended 
to  go  solely  into  the  controlled  modes  also  leaks  into  the  uncontrolled 
modes. 
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This  report  discusses  the  control  of  flexible  systems  described  by 
a  generalized  one-dimensional  wave  equation  which  relates  the  structure 
displacement  to  the  force  distribution  acting  on  the  structure.  Optimal 
control  involving  the  minimization  of  a  quadratic  performance  index  representing 
control  and  modal  energy  content  is  considered.  Typically  this  control 
formulation  leads  to  a  state  feedback  algorithm.  Since  the  state  components 
are  the  modes  of  vibration  which  are  not  directly  measurable,  a  means 
of  state  reconstruction  must  be  considered.  The  approach  taken  in  this 
report  is  a  deterministic  one  for  which  the  Luenberger  state  observer  is 
sufficient  for  state  estimation. 

II.  PROBLEM  FORMULATION 

The  structural  system  considered  in  this  analysis  is  expressed  in 
the  following  distributed  parameter  form 

~2 

— y(x,t)  +  A  u(x,t)  =  F(x, t)  (1) 

at 

where  u(x,t)  is  the  displacement  at  point  x  and  time  t,  A  is  a  linear 
self-adjoint  differential  operator  and  F(x,t)  is  an  externally  applied 
force  on  the  system  which  can  represent  the  control  force  or  an  external 
disturbance  or  both.  If,  for  instance,  the  control  is  provided  by  M  point 
actuators  located  at  x^,  1=1,2, -*-,M  and  there  is  no  other  external  disturbance, 
then 

M 

F (x , t )  =  l  <$(x-x,)q  (t)  (2) 

i=l 

where  5  is  the  Dirac  delta  function  and  q^(t)  is  the  force  applied  at 
System  Modal  Equations 

The  system  equations  of  motion  expressed  by  equation  (1)  and 
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accompanying  boundary  conditions  can  be  transformed  to  the  modal  representation 
by  use  of  the  expansion  theorem  [4] 


u(x,t)  =  l  <j>  (x)u  (t) 
r=l  r  r 


ur  (t) 


' 

u(x,t)<|ir  (x)dx,  r  =  1 , 2 ,  •  •  • 


(3) 

(4) 


where  u^ft)  are  modal  amplitudes  and  4>  (x)  are  the  eigenfunctions  of  the 
operator  A  appearing  in  equation  (1).  The  eigenfunctions  are  determined 
from  the  equation 

A$r(x)  =  Xr$r(x)  ,  r  =  1,2,  ••  •  (5) 

where  Xr  are  the  eigenvalues  of  the  operator  A.  Similarly,  the  control 
force  F(x,t)  can  be  represented  in  modal  form  as 

CO 

F(x,t>  =  l  ♦  <x)f  ft)  (6) 

r-1  r  r 

f r Ct )  =  |  F(x,t)<J>r (x)dx  (7) 


where  f^Ct)  are  the  modal  control  forces. 

The  modal  amplitudes  are  computed  from  the  following  equations  obtained 
by  substitution  of  the  results  of  equations  (3)  and  (6)  into  equation 
(1)  while  imposing  the  orthogonality  properties  of  the  eigenfunctions 


«r(t)  +  *rur(t)  =  f r (t )  ,  r  =  1,2, 


Initial  values  for  u^Ct)  are  obtained  from 


ur(0)  =  |  u(x,0)<f>r(x)dx 


ur(0) 


|~(x,0)<})r  (x)dx  . 


(8) 
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Modal  State  Equation 


A  larger  number  of  elastic  modes  may  be  needed  in  equation  (3)  for 
accurate  representation  of  the  displacement  of  a  large  flexible  structure. 
However,  the  concept  of  modal  control  is  to  restrict  active  control  to 
a  few  critical  modes  referred  to  as  controlled  modes.  The  remaining  modes 
which  are  needed  for  the  accurate  representation  of  the  structure  displacement 
are  referred  to  as  residual  modes.  If  we  choose  L  modes  for  descriptive 
purposes  but  only  control  N  we  have 

u(x,t)  =  uc<x,t)  +  UR  (x ,  t)  (9) 

where 

N 

uc(x,t)  =  51  <t>r(x)ur(t)  (10) 

r=l 

L 

u  (x,t)  =  l  <L(x)u  (t)  (11) 

r=N+l  r  r 

The  subscripts  u C lf  and  "R11  refer  to  controlled  and  residual  respectively. 

The  control  design  for  motion  suppression  based  on  only  N  modes,  when 
actually  L  modes  are  represented  in  the  physical  observations  and  are 
also  effected  when  the  control  is  activated,  can  lead  to  stability  problems 
referred  to  as  spillover.  This  phenomenon  is  discussed  in  some  detail 
in  references  [1,3,5], 

Equation  (8)  for  the  dynamics  of  the  controlled  modal  amplitudes 
can  be  expressed  in  state  space  form  as 

v(t)  =  Av(t)  4-  if(t)  (12) 

where 

f(t)  =  [f1(t)  f2(t)----fN(t)]T  (13) 
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v(t)  =  [u1(t)'*'uM<t)fi1(t)*"dH(t)]T 


B  = 


A 


where 


I 


N 


A 


1  0 
0  1 


0  0 

_ 


0 

0 


(NxN) 


1 


X1  0  •  •  •  0 

0  V  ’  *  0 


0  0 


X 


N 


(14) 

(15) 


(16) 


(17) 


(18) 


Control  Force  Description 

In  this  analysis  it  was  assumed  that  the  structure  control  could  be 
effected  by  M  point  actuators  placed  at  discrete  points  xi  along  the 
structure.  An  expression  for  this  force  was  given  in  equation  (2) 
and  is  repeated  here: 


M 

F(x,t)  =  5!  5(x-x  )q  (t)  (19) 

i=l 
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The  modal  force  is  computed  from  equation  (7)  as 


f  (t) 
r 


F(x,  t)<|)  (x)dx 


M  r 

=  I  ^(t)* 

i-lJ 


(tJdCx-x^dx 


M 


I  *  Cx±)q.(t) 
1-1 


(20) 


Substitution  of  f  (t)  as  given  by  equation  (20)  into  equation  (12)  results 


in  the  state  equation 

v(t)  *  A  v(t)  +  B  q(t) 


where 


q (t)  =  [qx(t)  q£(t)  qM(t)] 


(21) 


(22) 


and 


B  = 


<J>1(x1)'*^1(xM) 


(23) 


VX1)*‘‘VXM) 

Measurements 

For  purposes  of  feedback  control  it  is  asssumed  that  displacements 


y^(t)  -  u(Zj,t)  ,  j  =  1.2* • ■ * ,p 


(24) 


are  measured  at  P  points  z_,  along  the  structure.  The  modal  representation 
for  equation  (24)  is 
N 


y< (t)  =  l  u  (t)4  (z .)  +  r 
J  k=l  J  J 


(t) 


(25) 
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yCt)  =  [y1(t)  y2(t) • • *yp(t) ] 


T 


(26) 


then  we  have 


y (t)  =  C  v(t)  +  R(t) 


(27) 


where 


^  (z^)  •  •  ,<f>N(z1)  0*  *  *0 


C  = 


(28) 


vyvv  0-,-° 


and 


R(t)  =  [r1(t)  r2(t)  •••  rp(t)] 


T 


(29) 


III.  THE  CONTROL  PROBLEM 

In  this  study,  we  address  the  problem  of  vibration  suppression  in 
distributed  systems  governed  approximately  by  the  modal  equations  (21)  and 
(27).  We  assume  that  spillover  effects  are  negligible  and  set  R(t)=0  in  (27) 
The  active  control  of  large  space  structures  serves  as  our  primary  motivation 
In  these  structures,  typically  large  essential  masses  dictated  by  mission 
requirements  are  connected  by  low  mass,  flexible  trusses.  Normal  operations 
aboard  these  structures  naturally  introduce  vibrations  in  the  supporting 
trusses  which  must  be  damped  effectively.  Mathematically,  two  primary 
sources  of  vibrations  can  be  identified: 


(i)  Changes  In  initial  conditions  in  the  dynamics  (21)  (because 

of  slewing  or  fine  pointing  effects,  sudden  disturbances,  etc.); 


(ii)  Fixed,  known  disturbances  (because  of  onboard  manufacturing 


machinery,  pumps  and  motors,  etc.). 
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Each  of  these  sources  of  vibration  will  be  discussed  in  detail  below. 

An  initial  configuration  for  the  structure  governed  by  (1)  is  determined 
once  the  functions 

u(0,x)  and  |^<0 , x)  (30) 

are  specified,  since  these  determine  initial  conditions 


ur (0) ,  ur(0) ,  r  =  1,2, •  •  *  (31) 

for  the  modal  equations  (8).  Conversely,  initial  conditions  (31)  determine 
functions  (30).  Retaining  N  modes  to  approximate  (1),  as  in  (21),  appropriate 
system  initial  conditions  are  given  by  specifying  a  vector 


vQ  =  fv^O)  v2(0)  v2N(°)]T*  (32) 

The  suppression  of  vibrations  in  system  (21)  due  to  changes  in  initial 
conditions  (32),  as  described  in  (i)  above,  may  be  stated  precisely  as 
follows : 

Given  v(0)  =  v^,  find  control  q(t)  such  that  v(t)~K).  (33) 

Since  the  solution  u(x,t)  of  (1)  is  given  approximately  by  u  (x,t)  of 

c 

equation  (10),  if  v(t)-H),  Ug(x*t)-K)  also  for  each  x  and  vibrations  due 
to  the  initial  conditions 
N 

uc(x,0)  =  £  <t>r(x)vr(0)  (34) 

r=l 


and 


!^<x,°) 


N 


-  I 


r=l 


*rCx)vr+N(0) 


(35) 


are  suppressed. 

The  problem  in  (33)  above  is  customarily  viewed  as  a  linear  regulator 
problem  in  optimal  control  theory.  The  criterion  for  the  design  of  the 
control  algorithm  is  taken  to  be  the  minimization  of  a  quadratic 
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performance  index  which  is  a  function  of  energy  associated  with  the  N 
controlled  modes  and  the  control  effort  q(t).  More  precisely,  we  minimize 
over  all  controls  q  -  q(t) 

J(q)  ~  \  |  [vTQv  +  qTRq]dt  (36) 

1  J0 


subject  to 

v(t)  =  Av(t)  +  Bq (t)  (37) 

y(t)  *  Cv(t)  (38) 


where 


Q  - 


(39) 


A  as  in  (18),  and  R  is  a  positive  definite  weighting  matrix.  It  is  well- 
known  that  the  function  q*(t)  minimizing  (36)  is 


q* (t)  ■  -R  1BTSv(t) 


(40) 


where  S  is  an  NxN  symmetric,  non-negative  definite  solution  of  the  algebraic 


Riccati  equation 

SA  +  ATS  -  sbr“1bts  +  Q  =  0  (41) 

It  should  be  noted  that  in  practice,  in  order  to  implement  the  control 

A 

design  (40),  a  state  observer  v(t)  must  be  constructed  satisfying 


v(t)  -  Av(t)  +  Bq (t)  +  Gq [y (t)-y (t) ] 

(42) 

y(t)  =  Cv(t) 

(43) 

The  suboptimal  control 

q*(t)  =  -R  ^B^Sv(t) 

(44) 

is  then  used  to  actually  control  the  system  (37) . 
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We  now  consider  the  suppression  of  vibrations  due  to  fixed,  known 
periodic  disturbances,  as  described  in  (ii)  above.  If  M'  point  disturbances 
d^ (t) , • * • ,dM,  (t)  are  located  at  positions  x  =  x^,  i  =  1,2, along  the 
distributed  system  (1),  the  total  force  F(x,t)  acting  on  the  controlled 
system  is 

F  (x » t )  =  Fc(x,t)  +  FD(x,t)  (45) 

where  F„(x,t)  and  F_(x,t)  are  the  control  and  disturbance  forces  respectively. 

Lr  U 

As  in  (19),  we  may  write 

M  M' 

F(x, t)  =  l  5 (x-x . )q . (t)  +  l  6(x-x')d,(t)  (46) 

i=l  i=l 

The  approximate  system  dynamics  (21)  then  becomes 

v(t)  =  Av  ( t )  +  Bq (t)  +  B'd(t)  (47) 

where 

d(t)  =  [dj.Ct)  d2(t)---dM((t)]T  (48) 


and 


B'  =  ^(x')  •••  ^(x^,) 


(49) 


1J  IT  H#' 

To  control  (47),  we  make  use  of  a  very  simple  idea, 
a  control  gain  matrix 


Namely,  determine 


K  = 


11 


V 


kMl  •  •  •  •  ' 


(50) 
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such  that 


* 


q(t)  =  t q x (t )  q2(t)  ••**  ^ C t ) ] T,  (51) 

with 

qi(t)  “  Kid(t)  1 "  (52) 

suppresses  the  effects  of  vibrations  in  (47)  caused  by  d(t).  Ki  is  the 
ith  row  of  K.  (In  application,  d(t)  would  have  to  be  fed  back  with  a  slight 
time  delay.)  More  precisely,  the  vibration  suppression  problem  of  (ii) 
above  may  be  stated  as  follows! 

T 

Given  known  disturbances  d(t)  =  [d^(t) • • *d^, (t) ]  , 
find  K  such  that  when  (47)  is  controlled  by  q(t)  as 

(53) 

constructed  in  (51)  and  (52) ,  the  steady  state 
magnitudes  |v  (t)|  are  as  small  as  possible. 

To  minimize  the  effects  of  d(t),  one  should  attempt  to  minimize  the 
effects  of  the  vector  of  generalized  forces 


BKd(t)  +  B'd(t)  =  [BK+B']d(t)  (54) 

in  (47).  To  do  this,  we  attempt  to  minimize  the  magnitudes  of  the  elements 
in  the  matrix 

Bk  +  B'  <55> 


by  minimizing  the  quadratic  function 


F(K) 


1  2N  T 

I  l  W  (B  K+B')  (B  K+B')T. 
i=N+l 


(56) 


The  minimization  in  (56)  is  over  the  elements  k_  of  K  and  the  weights 
w^^  are  to  be  suitably  chosen.  F(K)  may  be  minimized  numerically,  or, 
the  linear  equations  that  result  when  its  gradient  is  equated  to  zero 
may  be  solved. 
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As  a  special  case,  to  reduce  computations  in  the  minimization  of 
F(K),  one  could  seek  a  vector 

K  “  [k1k2-.-kM,l  <57) 

that  minimizes  the  effects  of  d(t)  through  the  control  q(t)  with 

qi(t)  =  Kd(t)  i  =  1,2, *  * • ,M-  (58) 

In  this  simpler  case,  we  need  only  find  m'  constants  k^,  k2>  k  ,  which 

minimize  (56),  rather  than  the  MM'  constants  for  K  given  by  (50). 

IV.  A  NUMERICAL  EXAMPLE 

As  a  simple  illustrative  example  of  vibration  suppression  in  distributed 
parameter  systems,  we  consider  the  vibrational  control  of  the  cantilevered 
beam  pictured  in  Fig.  1.  The  partial  differential  equation  governing 
undamped  displacements  u(x,t)  of  points  x  (along  the  axis  of  the  beam) 
at  times  t  is 


^<x,t)  (x,t)  =  ^F(x,t) 

*  Z  m  _  4  m 

3t  3x 


with  boundary  conditions 


(59) 


2  3 

u (0 , t )  =  ~  (0,t)  =  ^  (L.t)  -  (L, t)  =  0 

3X  3x^  3xJ 


(60) 


Here  m  =  mass,  E  -  modulus  of  elasticity,  I  “  cross-sectional  area,  and 
L  =  length.  For  simplicity,  we  will  assume  that  m=E=I«L=l,  a 


single  point  accuator  is  located  at  x  -  — ,  and  (possibly)  a  single  disturbance 


acts  at  x  =  -j- 

Simulations  of  the  suppression  of  vibrations  in  the  beam  due  to  changes 
in  initial  conditions  and  known  external  disturbances  (as  discussed  in 
the  preceding  section)  were  performed.  Programs  were  coded  in  FORTRAN 
77  and  run  on  a  VAX  11-780  computer-  Time  responses  were  obtained  using 
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the  IBM  360  SSP  routine  DRKGS,  a  fourth  order,  variable  step  size  Runge- 
Kutta  integration  technique.  Algebraic  Riccati  equation  solutions  were 
obtained  using  a  collection  of  Kleinman  algorithms.  The  results  of  these 
simulations  will  now  be  discussed. 

Transforming  (59)  and  (60)  into  the  state  space,  retaining  two  modes 
for  analysis  (N=2),  and  introducing  internal  damping  via  the  damping  ratio 
C,  equations  (21)  for  the  cantilevered  beam  become 


r  0  0  1  0 

r0 

0  0  0  1 

0 

v(t)  = 

2 

0  -2£u)^  0 

v(t)  + 

<f>2%) 

0  -tiii?  0  -2?(ij- 

2  2J 

2  (3  )j 

q(t) 


(61) 


Ignoring  spillover  effects,  we  have  the  approximate  equality 
u(x,t)  ;  iji1(x)v1(t)  +  $2(x)v2(t) 


(62) 


where  in  general,  the  mode  shapes  ^(x)  for  (59)  are  given  explicitly  as 

a 

<J>  (x)  =  A  (sin  a  -sinh  a  )(sin  —  x-sinh  —  x) 


a  a 

r  ■£ 

+  (cos  a  +cosh  a  )  (cos  - —  x-cosh  —  x) 
r  r  L  L 


(63) 


for  explicitly  computable  constants  (see  [4]).  The  A^  are  normalizing 
constants.  The  associated  natural  vibrational  frequencies  satisfy 


a 


'El 


id  =  =  7^-  /  - —  . 

r  r  L  J  m 


EX  3 

(Note  that  A  = - r  in  equation  (5).) 

m  3x* 


(64) 


The  mode  shapes  <^(x)  and  ^(x)  are  represented  schematically  in  Fig. 

2.  For  all  simulations,  the  beam  is  assumed  configured  as  in  Fig.  1  with, 
as  previously  stated,  m=E=I=L=l.  We  also  set  5  =  .05. 
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We  first  consider  problem  (33),  the  suppression  of  vibrations  due 
to  changes  in  initial  conditions  (Recall,  d(t)=0).  An  effective  manner 
in  which  to  simulate  abrupt  changes  in  a  system's  initial  conditions 
is  to  excite  it  with  an  impulse  disturbance.  For  the  beam  of  Fig.  1, 
we  consider  excitation  forces  of  the  form 

F(x, t)  =  a  S(x-l)<$(t).  (65) 

Such  a  force  corresponds  to  an  impulse  of  magnitude  a  applied  to  the 
free  end  of  the  beam  at  t  =  0.  With  respect  to  the  state  space  coordinates 
v,  it  is  easy  to  check  that  excitation  with  (65)  is  equivalent  to  specifying 
initial  conditions 

VjCO)  =  v2(0)  =  0  ,  v^(0)  =  Ouf^Cl),  v2  =  a02(l)  .  (66) 

The  uncontrolled  time  response  of  (61)  excited  by  (65)  with  a  =  ,002  is 

given  in  Fig.  3.  Here  tip  displacements  are  plotted  against  elapsed  time. 

When  the  optimal  control  (40)  is  applied  to  suppress  the  vibrations  (here 

R  =  [1])  assuming  that  the  entire  state  v(t)  is  available  to  the  controller, 

the  vibrations  are  quickly  arrested  as  revealed  by  Fig.  4.  Assuming  that 

only  observations  of  v^(t)  and  Vj(t)  are  available  (C  =  [1100]  in  (38)), 

a  state  observer  v(t)  as  in  (42)  and  (43)  was  constructed.  Letting 
-1  T 

K*  =  R  B  S,  Gq  in  (42)  was  determined  so  that  the  eigenvalues  of 
A  -  GqC  were  positioned  at  fifty  times  the  real  part  of  that  eigenvalue 
of  A  -  BK*  with  smallest  real  part  (in  absolute  value).  The  response  of 
(61)  when  excited  by  (65)  and  controlled  through  (44)  is  given  in  Fig.  5. 

It  is  assumed  that  the  observer  initial  conditions  are 

^(0)  =  v2(0)  =  v3(0)  =  O4(0)  =  0 

so  that  incomplete  knowledge  of  v(t)  is  present  at  t  =0.  Some  degradation  in 
performance  is  obvious,  but  v(t)  effectively  converges  to  v(t)  after  t  =  1.5. 
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Consider  now  problem  (53)  in  which  an  external  disturbance  d(t)  is 
present  (refer  to  Fig.  1).  Let 

d(t)  =  g  sin  9t  (68) 

where  =  3.516  <  9  <  22.034  =  •  Tip  displacements  when  system  (61) 

is  driven  by  (68)  with  g  =  .03  from  zero  initial  conditions  (v(0)  =  0) 
are  pictured  in  Fig.  6. 

Beam  control  is  provided  by  q(t)  as  in  (58)  (M  =  1).  In  this  example, 

K  in  (57)  is  simply  a  scalar  which  is  found  by  minimizing  (56).  The  minimizing 
K  is  easily  found  to  be 

K  =-.717.  (69) 

Fig.  6  displays  the  tip  time  response  when  (61)  is  controlled  through 
(58)  and  (69).  Obviously,  this  simple  approach  to  vibration  suppression 
is  very  effective. 

Figures  7  and  8  show  beam  configurations  at  various  instants  (time 

2 

slices).  It  is  interesting  to  note  that  when  (68)  is  applied  at  x  =  -j, 
even  though  this  force  is  directed  in  an  upward  direction  initially,  the 
tip  is  thrust  downward.  This  clearly  displays  the  true  flexibility  of 
the  simple  model  chosen  for  analysis. 

V.  EXTENSIONS  AND  FUTURE  WORK 

The  analyses  presented  in  the  preceding  sections  can  be  extended  naturally 
in  numerous  ways,  and  the  authors  intend  to  pursue  this  work  in  the  future. 

In  particular,  attention  should  be  given  to  the  following  further  investigations: 

(1)  Increase  the  number  of  actuators,  sensors,  and  external  disturbances 
and  evaluate  the  consequent  effects  upon  system  performance; 

(2)  Compute  actuator  power  requirements  and  perform  total  control 
energy  sensitivity  analyses  with  respect  to  the  placement  of 
the  components  in  (1); 
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(3)  Consider  more  general  disturbances  d(t),  which  are  not  clearly 
identified,  via  spectral  analysis; 

(4)  Consider  the  effects  of  spillover  on  performance; 

(5)  Couple  the  ideas  of  (l)-(4)  with  FDL  finite  element  analysis 
software  to  study  realistic,  possibly  high  dimensional  models 
of  large  space  structures* 

Finally,  there  has  been  considerable  interest  recently  (e.g.,  see 

[6])  in  combined  structural  and  control  analysis*  Traditionally,  In 
aerospace  engineering,  control  system  design  has  been  completed  only 
after  the  controlled  vehicle  has  been  fully  designed.  Since  large  space 
structures  are  so  costly  to  deploy  and  since  in  practice  only  limited 
control  energy  is  available  for  maneuvering  and  vibration  suppression, 
the  idea  arises  that  perhaps  a  simultaneous  consideration  of  structural 
and  control  design  can  result  in  substantial  savings  in  cost,  complexity, 
and  mission  effectiveness*  With  respect  to  the  cantilevered  beam,  optimal 
controls  depend  directly  upon  the  structural  parameters  m,  E*  L,  and 
I.  The  variation  of  these  parameters  eventually  also  should  be  included 
in  the  analyses  of  ( 1 )  —  C 5 }  above,  and  ultimately,  again  using  finite 
element  methods,  the  authors  hope  to  extend  their  study  to  the  difficult 
questions  of  combined  structural  and  control  design  for  realistic  large 
space  structures. 
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Figure  4  -  Controlled  Tip  Displacement  -  Impulse  Excitation  and  Complete  Observations 
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Figure  5  -  Controlled  Tip  Displacement  -  Impulse  Excitation  and  Incomplete  Observations 


11(1, t) 


n 


23 


Figure  6  -  Uncontrolled  Tip  Displacement  -  Known  Disturbance  d(t)  =  .03  sin  9t 
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Figure  7  -  Controlled  Tip  Displacement  -  Known  Disturbance  d(t)  =  .03  sin  9t  {K  =  -.717) 
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Figure  8  -  Uncontrolled  Beam  Configurations  at  Selected  Times 
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Figure  9  -  Controlled  Beam  Configurations  at  Selected  Times  t  (K  =  -.717) 


